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O ■ Abstract 

Weak interaction contributions to hyperfine splitting and Lamb shift in light electronic and 
muonic atoms are calculated. We notice that correction to hyperfine splitting turns into zero for 
^* ; deuterium. Weak correction to the Lamb shift in hydrogen is additionally suppressed in comparison 

with other cases by a small factor (1 — 4 sin 2 Ow)- 
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I. INTRODUCTION 

Unexpected results of the Lamb shift measurement in muonic hydrogen at PSI [1] gave 
rise to the proton radius puzzle. It consists of the five sigma discrepancy in the value of the 
proton radius extracted from the muonic hydrogen experiment [l| on the one hand, and the 
values extracted from the electronic hydrogen (see review in the CODATA compilation 2j) 
and electron-proton S catte ri n g a on the other hand. A hn rst of theoretica! activity foUowed 
the PSI experiment. Old results on the Lamb shift in muonic hydrogen were recalculated 



and confirmed (see, e.g., |4| and references therein) ^proton structure and polarizability 
corrections were critically reevaluated and improved ( p, |6[ and references therein) , possible 
new physics explanations were explored ( [3, S] and references therein) . Despite all these 
efforts no resolution of the proton radius puzzle was found. It seems now that solution of 
this problem will require a lot of additional experimental and theoretical work, in particular 
precise measurements of different transition frequencies in muonic hydrogen and other light 
muonic atoms. An experimental program on measurement of transition frequencies in light 
muonic atoms is now in progress by the CREMA collaboration at PSI. First experimental 
data on hyperfine splitting (HFS) in muonic hydrogen is coming soon, and the results for 
muonic deuterium and helium ion are to follow [9]. In anticipation of these experimental 
data we calculate below weak interaction contributions to HFS and Lamb shift in light 



muonic atoms, generalizing old results for muonic hydrogen [10] (see also |11|). 

Effective low-energy field-theoretic weak interaction Hamiltonian due to neutral currents 

n, 

for the fundamental fermions has the form (see, e.g., |l2| ) 



Hz = ^7fl Sx ^ r(T3 ~ sin2 BwQ) ^) ' (1) 

where 9w is the Weinberg angle, Q is the charge operator in terms of proton charge, T3 = 
Ts(l — 7 5 )/2, T3 is the weak isospin, and summation goes over all species of fermions. 
Each current in this local four-fermion Hamiltonian contains a vector and an axial part. 
For nucleons axial parts are renormalized by strong interactions and should be multiplied 



by Qa — 1-27 (see, e.g., (I2j). Specializing for the case of lepton-nucleon interaction the 
Hamiltonian in Eq. ([1]) reduces to (an extra factor two arises because all fields enter each 
factor in Eq. (JT])) 



H z = -^= d 3 x [VwV^ - ^ (1-4 sin 2 Bw) ipi] 

2V2 J (2) 

X [gAi'nl^Ai - 4>nJ^n ~ 9 'A$ 'pi 'pY H 'p + 4>plv ( l ~ 4 sill 2 %) , 

where ^ is the lepton (electron or muon) field, and ip p and ip n are the proton and neutron 
fields, respectively. This Hamiltonian generates all weak interaction contributions considered 
below. 



II. WEAK INTERACTION CONTRIBUTIONS TO HYPERFINE SPLITTING 

The leading weak interaction contribution to HFS arises from interaction of axial currents 
in Eq. ($2§. In the leading nonrelativistic a ppr oximation only spatial components of axial 
neutral currents give nonzero contributions [lOj, and the Hamiltonian simplifies 



H z / d 3 x $a^%) $„7M7fyn - i> P l»l%) 



2y/2 
9aGf 



d 3 X (iprfj 5 ^) (^VVn ~ Ippf^lpp) 



(3) 



2^2 

For a nucleus with Z protons and A — Z neutrons this field-theoretic Hamiltonian in the 
nonrelativistic limit reduces to the quantum mechanical Hamiltonian 

Hz = 9 £=-*>- (>>,->>, KV). (1) 



V P n / 



2s/2 

Matrix elements of this operator give the leading weak interaction contributions to hyperfine 
splitting that is nonzero only in S states. The only remaining task is to calculate expectation 
value of the scalar product taking into account the nuclear wave function. We consider below 
in parallel light electronic and muonic atoms and ions, but numerical results are provided 
only for muonic systems. 



A. Hydrogen 

In the case of muonic (electronic) hydrogen there is only one term in the nuclear factor 
in Eq. 01]), and we immediately obtain the leading weak interaction contribution to HFS 



splitting in the nS'-state 



lOj in the form 



AE z (nS) = 9 -0^\MO)\ 3 (<re ■ <t p )\ F fZI (5) 

where i^ n (0) = a/ (Zam r ) 3 / (7m 3 ) is the Coulomb-Schrodinger wave function at the origin 
[Z = 1 for hydrogen), m r = mim p /(mi + m p ) is the reduced mass, J = <Ti/2 is the lepton 
spin operator, i" = <t p /2 is the proton (nucleus) spin operator, F = I + J is the total angular 
momentum. Obviously, (cr e • cr p )|^=J = 4 and 

AE z (nS) = ^ A G F {Zamf _ 
V2 

Numerically for n = 2 in muonic hydrogen the weak contribution is 

AE Z (2S) = 2.8 x 1(T 4 meV, (7) 

what is at least an order of magnitude smaller than the uncertainty in hyperfme splitting 
due to proton structure contributions [13) ]. 

To elucidate the magnitude of the weak interaction contribution let us compare it with 
the dominant Fermi contribution to hyperfme splitting in muonic hydrogen is (see, e.g., 

HQ) 



E F = V (ZQ)3m - - 182.44 meV, (8) 
3 vrt\m v 

where g p ~ 5.58 ... [2| is the proton ^-factor in nuclear magnetons. 

For excited states the dominant contribution to hyperfme splitting scales as 1/n 3 and for 
the state with an arbitrary principal quantum number n the ratio of the weak and dominant 
contributions to HFS in muonic hydrogen is 



n 3 AEz(nS) 3 gAGprn^mp 
E F 2^/2,71 g P ot 



1.2... x HT 5 . (9) 



B. Deuterium 



Deuteron is a spin one loosely bound system of two nonrelativistic nucleons that are 
predominantly described by the ^-state wave function. Respective spin wave function is 
symmetric and the matrix element of the spin-symmetric deuteron nuclear factor of the 
effective Hamiltonian in Eq. (jl]) in this approximation is equal zero, 



(a p - cr n ) = 0. (10) 

This conclusion remains valid even with account of the admixture of the D wave in the 
deuteron wave function, since the D wave spin function is also symmetric with respect to 

n 

spin variables (see, e.g., [16j). Hence, the weak interaction contribution to hyperfine splitting 
in electronic and muonic deuterium in the leading nonrelativistic approximation is zero. 

C. Tritium 

Triton is a spin one half (I — 1/2) system of one proton and two neutrons (Z = 1, A = 3). 
The third component of isosp in for the triton is minus one half (T 3 = 1/2—1/2—1/2 = —1/2). 
It is predominantly described by a product of the S wave coordinate wave function and 
a completely antisymmetric spin-isospin wave function. Obviously, in this approximation 
(cr p — er ni — cr n2 ) = 21. A more accurate analysis with account for other components of the 
triton wave function produces [17J 

(a p - tr ni - cr n2 ) =21 U- ^P s > - \pd \ = 2d, (11) 

where c ps 0.92. 

Further calculations go exactly like in the hydrogen case above and we obtain 



2cqaGf (Zam r 



\3 



AE z (nS) = (12) 



y/2 nn 3 

Numerically for n = 2 in muonic tritium the weak contribution is 



AE z (n = 2) = 3.1 x 10~ 4 meV. (13) 

Like in the hydrogen case we compare the weak contribution with the dominant Fermi 
energy in muonic tritium 

Ep = -g t -± '- r - = 239.919 . . . meV, 14 

3 mim p 

where g t = 5.957924896(76) [2] is the triton g-factor in nuclear magnetons. 



5 



For excited states the dominant contribution to hyperfine splitting scales as 1/n 3 , and for 
the state with an arbitrary principal quantum number n the ratio of the weak and dominant 
contributions to HFS in muonic tritium is 



n 3 AE z (nS) = 3 cg A G F m^m p _ ^ Q ^ 1Q _ 5 
E F 2V2it g t Oi 



D. Helium Ion 

Helion is a spin one half (J = 1/2) system of two protons and a neutron [Z = 2, A = 3). 
The third component of isospin for the helion is one half (T3 = 1/2 + 1/2 — 1/2 = 1/2). 
Like the triton the helion is predominantly described by a product of the S wave coordinate 
wave function and a completely antisymmetric spin-isospin wave function. Obviously, in 
this approximation (a Pl + cr P2 — cr n ) = —21. A more accurate analysis with account for 



other components of the helion wave function produces 17] 



(tr Pl - ct P2 - a n ) = -21 ^1 - -P S , - -P D j = -2c/. (16) 
Further calculations go exactly like in the hydrogen and tritium cases above and we obtain 

AE z (nS) = _ 2c 9AG F (Zamr)\ 
y/2 vrn 3 

Numerically for n = 2 in muonic helium the weak contribution is 



AE z (n = 2) = -2.5 . . . x 1(T 3 meV. (18) 

Like in the hydrogen and tritium cases we compare the weak contribution with the dom- 
inant Fermi energy in muonic helium 

E F = = -1370.8... meV, (19) 

3 m„m n 



where gn = —4.255250613 J2] is the helion (^-factor in nuclear magnetons. 

For excited states the dominant contribution to hyperfine splitting scales as 1/n 3 , and for 
the state with an arbitrary principal quantum number n the ratio of the weak and dominant 
contributions to HFS in muonic helium is 
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E. Helium Ions e 4 He + , iu. 4 He + 



Spin of a-particle is zero and there is no hyperfine structure in e A He + 1 [i A He + helium 
ions, and respectively no weak interaction contribution to hyperfine structure. 

III. LEADING WEAK INTERACTION CONTRIBUTION TO LAMB SHIFT 

The leading weak interaction contribution to the Lamb shift arises from interaction of vec- 
tor currents in Eq. ([2]). In the leading nonrelativistic approximation only time components 
give nonzero contributions [10!], and the interaction Hamiltonian simplifies 

G f 

H z ->■ -4= / d 3 x (1-4 sin 2 6 W ) fa) (V»n7^n ~ ^ v l^ (1-4 sin 2 9 W ) ip p ) 

G F f - , , (21) 

~* 2^h I d3x ^ _ 4 sin2 ^ ^ n7 °^ n ~ ^ p7 ° ( x _ 4 sin2 ® w ) ^pp) ■ 

For a nucleus with Z protons and A — Z neutrons this field-theoretic Hamiltonian in the 
nonrelativistic limit reduces to the quantum mechanical Hamiltonian 

H z = ^= (l-4sin 2 ^) (A-Z-Z(l-4sin 2 9 w ))5^(r). (22) 
2v 2 

Matrix element of this operator gives the leading weak interaction contributions to the Lamb 
shift that is nonzero only in S states. We obtain an explicit expression for the leading weak 
correction in all light atoms in the form 

AE z (nS) = (1 - 4 sin 2 6 W ) [(A - Z) - Z (l - 4 sin 2 9 W )} ^ Z f\ ( 23 ) 
2\/2 i — I 

For A = Z = 1 this result was obtained in 10|. It is interesting to notice that in muonic 



hydrogen due to A = Z = 1 the weak contribution to the Lamb shift is additionally sup- 
pressed by a small factor 1 — 4 sin 2 6y/ ~ 0.08. This suppression disappears for all other light 
muonic systems. 

Let us compare weak contribution to the Lamb shift with the dominant contribution. The 
principal contribution to the Lamb shift in light muonic atoms is generated by the diagram 



with the electron vacuum polarization insertion in the Coulomb photon, and was calculated 
long time ago [18| (see also reviews in [14j, |l5j, Il9| ) 



where 



A ^ = _8a^ g2)(/3)> (24) 



Q$(P) - f Pdp J™ d(f nl (£) e"^ (l + ^ V^i, (25) 
p\ _ /(n-Z - 1)! /2pV _ £ 2m /2p 



L n J^_ 1 (x) is the associated Laguerre polynomial 20] , and /3 = m e /(m r Za) 



For the experimentally relevant interval 2P — 2S we obtain 

AE(2P - 25) = A£ 21 - AE 20 = ^f^ (Q$(P) - Q$(/3)), ( 2 ?) 

and 

AE Z (L,2S) 3G F m 2 r Z (1- Asm 2 6 W ) [(A- Z)- Z (1- Asm 2 9 W )] 

^ p ~ 25) ~ W^fl-gg^)) ' 

This expression demonstrates once again that due to the condition A = Z = 1 the weak 
interaction contribution to the Lamb shift is additionally suppressed by an extra factor 
1 — 4 sin 2 6\y ~ 0.08 in comparison with the weak interaction contribution in other light 
muonic systems. For muonic hydrogen /3 ~ 0.7, Q%j(ft) — 0.056, Q^iifi) — 0.0037, and we 
obtain 

AE z (L,n = 2) q . , 

4^ r « -1-7 x 10~ 9 . 29 

AE(2P - 25) y ! 

We see that the weak correction to the Lamb shift in muonic hydrogen is orders of magnitude 
smaller than relative error of the Lamb shift measurement ^ is also much smaller than 
uncertainties of the proton structure corrections |||. 



IV. CONCLUSIONS 



We calculated the leading weak contributions to HFS and Lamb shift in light muonic 
atoms and ions. The leading correction to HFS in deuterium is zero because the deuteron 

8 



weak interaction Hamiltonian is antisymmetric with respect to nucleon spin variables while 
the deuteron spin wave function is symmetric. Corrections to Lamb shift in hydrogen are 
additionally suppressed by the small factor (1 — 4 sin 2 9 W ). This happens because all other 
nuclei contain neutrons that weakly interact with leptons without this suppression factor. 
In all cases weak corrections are much smaller than current experimental and theoretical 
errors. 
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